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The correlation between level velocities and eigenfunction 
intensities provides a new way of exploring phase space lo- 
calization in quantized non-integrable systems. It can also 
serve as a measure of deviations from ergodicity due to quan- 
tum effects for typical observables. This paper relies on two 
well known paradigms of quantum chaos, the bakers map and 
the standard map, to study correlations in simple, yet chaotic, 
dynamical systems. The behaviors are dominated by the pres- 
ence of several classical structures. These primarily include 
short periodic orbits and their homoclinic excursions. The 
dependences of the correlations deriving from perturbations 
allow for eigenfunction features violating ergodicity to be se- 
lectively highlighted. A semiclassical theory based on peri- 
odic orbit sums leads to certain classical correlations that are 
super-exponentially cut off beyond a logarithmic time scale. 
The theory is seen to be quite successful in reproducing many 
of the quantum localization features. 

PACS numbers: 05.45.-a, 03.65.Sq, 05.45.Mt 



I. INTRODUCTION 

For a bounded, classically chaotic system, ergodicity is 
defined with respect to the energy surface, the only avail- 
able invariant space of finite measure. In an extension de- 
veloped just over twenty years ago, the consequences of 
ergodicity for the eigenstates of a corresponding quantum 
system were conjectured to give rise to a locally, Gaussian 
random behavior [|,@]. Shortly thereafter, work ensued 
on defining the concept of eigenstate ergodicity within 
a more rigorous framework 0. Some of the paradoxes 
and peculiarities have been recently explored as well . 
One expression of eigenstate ergodicity is that a typi- 
cal eigenstate would fluctuate over the energy surface, 
but otherwise be featureless, in an appropriate pseudo- 
phase-space representation such as the Wigner transform 
representation Q. Any statistically significant deviation 
from ergodicity in individual eigenstates is termed phase 
space localization. 

It came as a surprise when Heller discovered eigen- 
states were "scarred" by short, unstable periodic or- 
bits |||7]]. A great deal of theoretical and numerical re- 
search followed [p[-p^[, and experiments also [ p^Jl5| . In 
fact, scarring is just one of the means by which phase 
space localization can exist in the eigenstates of such sys- 
tems. Another means would be localizing effects due to 
transport barriers such as cantori ]Tfj| , |T7| ] or broken sepa- 
ratrices Despite these studies and the semiclassical 



construction of an eigenstate |lfj| ] , the properties of indi- 
vidual eigenstates remain somewhat a mystery. 

Individual eigenfunctions may not be physically very 
relevant in many situations, especially those involving 
a high density of states. In this case, groups of states 
contribute towards localization in ways that may be un- 
derstood with available semiclassical theories. One sim- 
ple and important quantity where this could arise is the 
time average of an observable as this is a weighted sum 
of several (in principle all) eigenfunctions. In the Heisen- 
berg picture, where the operator is evolving in time, the 
expectation value of the observable could be measured 
with any state. Phase space localization features would 
be especially evident if this state were chosen to be a 
wave packet well localized in such spaces. 

One of us has already studied the correlation between 
level velocities and wavefunction intensities in connection 
with localization ]2(J . This can be directly connected to 
the issues raised above, and our treatment thus extends 
the previous work. The present paper is a companion to 
a study of similar problems in continuous Hamiltonian 
systems (as opposed to Hamiltonian maps, i.e. discrete 
time systems) and billiards [ plj . The methods used in 
these two papers complement each other and the results 
in the present paper are detailed as the systems studied 
are much simpler. 

In a system with a non-degenerate spectrum the time 
average of an observable A in state |a) is 



A(t)\ = J2\(^n)\ 2 ^n\A\^ n ) 



(1) 



where \ip n ) are the eigenstates of the Hamiltonian H. 
If the system depends upon a parameter A which varies 
continuously, an energy level "velocity", dE n /d\, for the 
n th level can be defined (velocity is a bit of a misnomer 
for it is actually just a slope - we are not evolving the 
system parameter in time). The level velocity-intensity 
correlation is identical to the time average above if we 
identify dH/dX with A, as from the so called Hellmann- 
Feynman theorem (for instance |22[): 



(2) 



It must be noted that while we call the operator aver- 
age Eq. (|l|) a "correlation" it is not the true correla- 
tion that is obtained by dividing out the rms values of 
the wavefunction intensities and the operator expectation 
value (as defined in [E0[). In other words we are going 
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to study the covariance rather than the correlation. This 
is followed in this paper for two reasons; first, dividing 
out these quantities does not retain the meaning of the 
time average of an observable and second, the root mean 
square of the wavefunction intensities which is essentially 
the inverse participation ratio in phase space is itself a 
fairly complex quantity reflecting on phase-space local- 
ization. 



II. CORRELATION FOR THE QUANTUM 
BAKERS MAP 

A. Generalities 

We first formulate in general the quantity to be stud- 
ied and the approach to its semiclassical evaluation. The 
classical dynamical systems are discrete maps on the 
dimensionless unit two-torus whose cyclical coordinates 
are denoted (q,p). The quantum kinematics is set in a 
space of dimension N ]2^-|25[| where this is related to 
the scaled Planck constant as N = 1/h, and the clas- 
sical limit is the large N limit. The quantum dynam- 
ics is specified by a unitary operator U (quantum map) 
that propagates states by one discrete time step. The 
quantum stationary states are the eigensolutions of this 
operator. The N eigenfunctions and eigenangles are de- 
noted by {\ipi) , (pi', i = 0, • • • ,N — 1}. The eigenval- 
ues lie on the unit circle and are members of the set 
{exp(-i<f>i)i = 0,---,N-l}. 

The central object of interest is: 



C A (a) 



'a\U- n AU n \a) 



(4) 



N-l 



c A (a) = \{*m\ 2 mm 



(3) 



The operator A is a Hermitian operator and either de- 
scribes an observable or a perturbation allowing one to 
follow the levels' motions continuously. The state |a) 
represents a wave packet on the two-torus that is well 
localized in the [q,p) coordinates We will be in- 

terested in the quantum effects over and above the clas- 
sical limit and we will require that the operator is trace- 
less. Otherwise we will need to subtract the uncorrelated 
product of the averages of the eigenfunctions (unity) and 
the trace of the operator. This immediately also implies 
that the correlation according to Random Matrix The- 
ory (RMT) [^7j is zero as well. The ensemble average 
of Ca{oi) will wash out random oscillations that are a 
characteristic of the Gaussian distributed eigenfunctions 
of the random matrices. Specific localization properties 
that we will discuss are then not part of the RMT models 
of quantized chaotic systems. In the framework of level 
velocities we are considering the situation where the av- 
erage level velocity is zero, i.e., there is no net drift of 
the levels. 

As we noted in the introduction the correlation C A {a) 
is simply the time average of the observable: 



where the large angular brackets denotes the time (n) av- 
erage. This requires that degeneracies do not exist and 
we assume that this is the case as we are primarily in- 
terested in quantized chaotic systems. A physically less 
transparent identity that is nevertheless useful in subse- 
quent evaluations is: 

C A {a) = ((a\U n \a) Tr (it/-") ) (5) 

This may be written more symmetrically as 



C A (a) = (Tr (\a)(a\U n ) 



Tr AU~ 



(6) 



Thus, the correlation is a sort of time evolved average 
correlation between the two operators A and |a)(a|. The 
semiclassical expressions for these are however different 
as complications arise from the classical limit of |a)(a| 
which would be varying over scales of H that govern the 
validity of the stationary phase approximations. How- 
ever, we may anticipate, based on the last form, that 
the semiclassical expression would be roughly the corre- 
lations of the classical limits of these two operators Eq] . 



B. Semiclassical Evaluation 

The bakers map is a very attractive system to study the 
quantities discussed in the introduction. The classical dy- 
namics is particularly simple (it is sometimes referred to 
as the "harmonic oscillator of chaos"). A simple quanti- 
zation is due to Balazs and Voros |25| (where a discussion 
of the classical dynamics may also found). As a model of 
quantum chaos it shows many generic features including 
the one central to this study namely scarring localization 
of eigenfunctions There are detailed semiclassical 

theories that have been verified substantially @,^9| 
We neglect certain anomalous features of the quantum 
bakers map 1 29 3^] that would eventually show up in the 
classical limit. This is reasonable in the range of scaled 
Planck constant values we have used in the following. 

We use the second time averaged expression, Eq. (|^), 
for the correlation. We do not repeat here details of 
the quantization of the bakers map or the semiclassical 
theories of this operator except note that we use the anti- 
periodic boundary conditions as stipulated by Saraceno 
]l2] | in order to retain fully the classical symmetries. 

The semiclassical theory of the bakers map deals with 
the powers of the propagator. The trace of U n , the time 
n propagator, has been written in the canonical form 
of a sum over classical hyperbolic periodic orbits with 
the phases being actions and the amplitudes relating to 
the linear stability of the orbits. The complications with 
Maslov phases is absent here |2^,^]. Also, the semiclas- 
sical expressions have been derived for matrix elements 
of the time n propagator in the wave packets basis p6[. 
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The time domain dominates the study of the quantum 
maps, the Fourier transform to the spectrum being done 
exactly. Our approach to the correlation is then natu- 
rally built in the time domain. The situation is different 
in the case of Hamiltonian time independent flows where 
the energy domain is very useful. 

We use the semiclassical expression for the propagator 
diagonal matrix elements derived in ]2q] : 



1*7" I 



E 



exp(iS-y/K) 



E 



exp 



. ^cosh(An) ( 

iSqSp 



cosh(nA) — 1 
2 cosh(rtA)?i 



x(<v + V) 



tanh(nA) 



(7) 



Here 7 labels periodic orbits of period n including repe- 
titions. The Lyapunov exponent is A which is ln(2) for 
the usual bakers map (corresponding to the (1/2,1/2) 
partition and Bernoulli process). Also Ti = h/(2w) — 
l/(2wN), 5q 

— 1j ~ Qa and a similar relation for p. The 
position of j-th periodic point on the periodic orbit 7 
is (qj,pj). The centroids of the wave packets, assumed 
circular Gaussians, are (q a ,p a ). The choice of type of 
wave packets is not crucial for the features we seek. We 
note that the simplicity of this expression for the prop- 
agator derives from the simplicity of the classical bak- 
ers map, especially the fact that the stable and unstable 
manifolds are everywhere aligned with the (q,p) axes. 
That Eq. (0) happens to be a periodic orbit sum differs 
from the similar treatment for billiards as found in f32| | 
where such sums are treated as homoclinic orbit sums. 
Note however, that the local linearity of the bakers map 
renders the two approaches (periodic orbit, homoclinic 
orbit) equivalent. 

A generalization of the trace formula for the propa- 
gator is given below that is easily derived by the usual 
procedure employed for the propagator itself |3FJ. Such 
a formula was derived in |53j for the case of Hamiltonian 
flows in the energy domain. We make the simplifying 
assumption that the operator A is diagonal in the posi- 
tion representation (we could treat the case of A being 
diagonal in momentum alone as well). This avoids the 
problem of a Weyl-Wigner association of operators to 
functions on the torus. The quantum operator A under 
this simplifying assumption has an obvious classical limit 
and associated function which is denoted by A(q). The 
other major assumption used in deriving the formula be- 
low is that it does not vary on scales comparable to or 
smaller than Ti. 

Thus we derive: 



Tr (if/"") ~ 



exp(— iS^/K) 



2sinh(nA/2) 



E^fe). 



(8) 



The index j again labels points along the periodic orbit 
7. The sum over the periodic orbit is the analogue of the 
integral of the Weyl transform over a primitive periodic 
orbit in the Hamiltonian flow case [B3| . The special case 



A = I the identity corresponds to the usual trace formula 
p9|^0[ . Note that we have written the sums above as 
being over periodic orbits, while the trace formulas have 
been often written as sums over fixed points. 

The first step is to multiply the two semiclassical peri- 
odic orbit sums in Eq. (Q) and in Eq. (||). Since there is 
a time average, n is assumed large enough, but not too 
large (so that these expansions retain some accuracy). 
All hyperbolic functions are approximated by their dom- 
inant exponential dependences. The diagonal approxi- 
mation and the uniformity principle p4J is used as well. 

C A (a) = {J2V2cM-n\)J2 ( E F fe,^) 



'A(T qj ,T Pj 




(9) 



Here we have taken a more general dependence for A 
(including the possibility of momentum dependence). T 
represents elements of the symmetry group of the sys- 
tem including time-reversal symmetry and including, of 
course, unity. These symmetries imply in general, though 
not as a rule, distinct (for T ^ I) orbits with identical 
actions. One assumes that the overwhelming number of 
action degeneracies are due to such symmetries. 
The function F is the approximated Gaussian: 

iSqSp 



exp 



1 

2h 



(Sq' + Sp 2 ) 



(10) 



Using A = ln(2) and the fact that there are approximately 
2 n /n orbits of period n, one finds 



C A (a) = V2~Y. 



M 

E 

T l = -M 



C T (l) 



where C(l) is a classical Z-step correlation: 
1 - 

Ct{1) = -J2F( qj , Pj )A(Tq j+l ,Tp j+l ). 



(11) 



(12) 



3=1 



The time average is taken over a typical orbit. We aban- 
don any specific periodic orbit and appeal to ergodicity, 
taking n and also M as practically infinite. This is with 
the assumption that such correlations will decay with 
time I. In fact, below we calculate such correlations ex- 
plicitly and display the decay. Note that Ct{1) 7^ Cr(— Z) 
in general. Although these are classical correlations, in 
the sense that qj,Pj represent a classical orbit, h appears 
as a parameter in them through F. Further, using the 
ergodic principle we can replace time averages in Ct(1) 
by appropriate phase space averages: 



dq dpF(q,p)A(Tf l (q,p),Tg l (q,p)) (13) 



C T (l) 



where we have used the fact that the total phase space 
volume (area) is unity, and f l (q,p) = qi, g l {q,p) = Pi are 
the classical Z-step integrated mappings. 
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C. Special case and verifications 

We first consider the case that A = A (T p + T$)/2, 

where T p is the unitary single-step momentum transla- 
tion operator that is diagonal in the position representa- 
tion and Ao is a constant real number. This implies that 
the associated function is A{q) = A cos{2irq). Below 
we consider Aq = 1 as the strength of the perturbation. 
The elements of T, apart from the identity (J), are time- 
reversal (TR) symmetry and parity (P). Time reversal 
in the bakers map is (T(q) — p, T{p) — q) followed by 
backward iteration, while parity is the transformation 
(T(q) = l-q,T(p) = l-p). 

We begin with the evaluation of the forward correlation 
(I > 0) corresponding to T = I. 



dqdpF(q,p) cos(2ir2 l q). 



This follows from the equality: 

f\q) = 2<< 7 (modl) 



(14) 



(15) 



for the bakers map. The limits of the integrals can be 
extended to the entire plane as long as the centroid of 
the weighting factor (q a ,p a ) is far enough away from the 
edges of the unit phase space square that the Gaussian 
tails are small there. The integral is elementary, and 
using h = 1/N one gets: 



C^l > 0) 



1 



V2N 



exp( 



-2 2l ir/(2N))cos(2Tr2 l q a ). (16) 



This explicit expression shows the super-exponential de- 
crease with time I in the correlation coefficients. It 
is interesting to note that the logarithmic time scale 
which sets an important quantum-classical correspon- 
dence scale of divergence for chaotic systems, here r = 
1/A Ln(l/27r7i) = ln(A)/ln(2), enters the correlation de- 
cay. In fact, the correlations are significant to precisely 
half the log-time. We anticipate this feature to hold in 
general, including autonomous Hamiltonian systems. 

Since g~ l (p) = 2'p(mod 1), for (I > 0), the time- 
reversed backward correlation (I < 0) is 



C T r{1 < 0) 



1 



exp(-2- 2i 7r/(2A0) cos(2tt2-V) 



V2N 



(17) 



which also decays super-exponentially and is responsible 
for the (q <-> p) symmetry in the final correlation. 

Next we turn to the other, apparently more curious, 
correlations: the backward identity correlations and the 
forward time-reversed one. As an example of a backward 
(I < 0) identity correlation consider I = — 1: 



q/2 for p< 1/2 

(q + l)/2 forp > 1/2 



Therefore 
Ci(-l) 



1 ,1/2 

dq / dpF{q 1 p) cos(7rg) 
Jo 

dq dpF(q,p)cos(irq) 

JO Jl/2 



(18) 



In fact, since cos(7rg) vanishes at 1/2, there is no dis- 
continuity in the full integral, but it is more difficult to 
evaluate (and to approximate). If one were to take the 
upper limits of the p integrals to be infinity, there would 
be errors at p = 1/2. However, this is not terribly dam- 
aging, and tolerating a small discontinuity at this point 
due to this approximation leads to: 



Ci(-1) 



± 



1 



V2N 



exp(— 7r/(8JV)) cos(nq 



(19) 



the sign depending on if p a < 1/2 or if p a > 1/2 respec- 
tively. The time-reversed, forward correlation, Ctr(X), 
is the same as this except for interchanging the roles of 
q a and p a . 

The generalization of this to higher times is (take I > 
below) : 



Ci(-l) = Y, I d< lf ( + >/ d P F(q,p)cos(27r(q + V)/2 l ) 

V=Q Jo Ju/2' 

(20) 



where v represents a partition of the bakers map at time 
I, and V results from the bit-reversal of the binary expan- 
sion of v. The momentum gets exponentially partitioned 
with time, and it precludes going beyond the log-time 
here as well (like the forward correlation) , although there 
is apparently no super-exponential decrease here. Indeed 
if we evaluate the above after neglecting finite limits in 
each of the p integrals above, so that we would have 2 l 
discontinuities at time I, we get: 



Ci(-l) 



1 



V2N 



exp(-7r/(2( 2i+1 ) AO) cos(27r(g Q + 77) /2 l ) 

(21) 



depending on if p a lies in the interval (v/2 l , [y + l)/2 l ). 
So that for I large and N fixed, the exponential goes to 
unity; effectively, for large N and any I, the exponential 
can be replaced by unity. Even the g Q -dependent part 
of the argument in the function (cos) itself is tending to 
vanish, so that the integral seems to give the area of the 
Gaussian (h). The approximation of putting all p lim- 
its to infinity makes sense only if the Gaussian state is 
well within a zone of the partition and this is necessarily 
violated at half the log-time. So the approximate expres- 
sion of Eq. ( pl|) breaks down beyond t/2. This lack of a 
super-exponential cutoff as seen with the previous corre- 
lations considered is due to two special conditions. First, 
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the argument of the cosine has no p-dependence. Sec- 
ond, all the stable manifolds are perfectly parallel to the 
p axis. We would recover super-exponential decay in all 
the correlations if the operator, A, being considered was a 
constant function along neither the stable nor the unsta- 
ble manifold. In this sense, we have chosen a maximally 
difficult operator with which to test the scmiclassical the- 
ory, though it simplifies the quantum calculations. 

As before, Ctr{1){<1,p) = Ci(—l)(p,q). Parity symme- 
try is benign and leads to an overall multiplication by a 
factor of 2. Thus, the final scmiclassical expression for 
the full correlation for the quantum bakers map is: 



C A {a) 



N 



Ti 



exp(-2 2/ 7r/(2A0) cos(27t2 z <7 q ) 



1=0 



T-2 



2-1 



+ 5>xp(-^/(2 2J+1 A0) ]T (cos(27r(g Q +V)/2 l ) 

1 = 1 !/=0 

xe(p Q -^/2')e((^ + i)/2 i -p Q ))] 

+ {q a ^p a )- (22) 

where T\ can be infinite but it is sufficient to stop just 
beyond half the log-time. As just discussed, T 2 , is more 
problematic here, and we do not have an expression to 
use beyond r/2. 9 is the Heavyside step function that 
is zero if the argument is negative and unity otherwise. 
The correlation is of the order 1/TV or h. If one were to 
divide by the number of states in Eq. (^) so that it is 
a true average, this quantity would decrease as 1/TV 2 or 

h 2 . 



sum (over T% terms). Figure 1(b) compares the semi- 
classical formula to the exact quantum calculation. We 
have taken eight "forward" correlations (excluding zero) , 
i.e. Ti = 8, while we have only taken two "backward" 
correlations, i.e. Ti = 2. This is because it appears 
that the approximations that go into the latter expres- 
sions lead to non-uniformly converging quantities and it 
works better to stop at a earlier point in the series. The 
(artificial) discontinuities at 1/2 and 1/4 are seen promi- 
nently in the semiclassical results. Otherwise, it turns 
out that the semiclassical approximation captures many 
fine-scale features of the correlations, some of which will 
be discussed below. Figures (2a, b) are for specific one 
dimensional sections of the same quantities. The agree- 
ment is very good. 





< 




FIG. 2. Sections of the correlation for N = 100. (aj 
Pa = -33 section, (b) p a = .72 section. The points are the 
quantum calculation while the solid lines are semiclassical 
evaluations. 



FIG. 1. (a) The absolute value of the quantum correlation 
with the cos(27rg) operator for TV = 100 whereas (b) corre- 
sponds to a semiclassical evaluation of the same. 



For the case of TV — 100, we compare 



in Fig. (1) the 
with the final 



full quantum correlation given by Eq. 
semiclassical evaluation given by Eq. (E2 ) . The absolute 
value of the correlation function is contoured and super- 
posed on a grey scale. Figure (la) shows the quantum 
calculation for the full phase space. In other words, the 
intensity (value) of each point, (q,p), on the plot repre- 
sents the CA(o;)-calculation for a wave packet centered 
at (q a = q,p a = p)- The first sum in Eq. ( p2[ ) (over T\ 
terms) is a smooth function, and it also displays an ad- 
ditional symmetry about 1/2 in both canonical variables 
separately. This extra symmetry is broken by the second 



D. Classical features in the correlation 

A strong (positive) correlation is indicated at the clas- 
sical fixed points (0, 0) and (1, 1), with the rest of the sig- 
nificant correlations being negative. They are dominated 
by several classical structures as illustrated in Fig. (3). 
Here the TV value used is 200, and superposed on the sig- 
nificant contour features are the following classical orbits: 

i) the period-2 orbit at (1/3, 2/3), (2/3, 1/3) is by far 
the most prominent structure. This is shown in 
Fig. (3a) by two circular dots. Also, we can look 
at these structures closely through 1-d slices. In 
Fig. (2a), the correlation is seen to be large and 
negative at (q a — .66). The period-2 structure is 
dominating the landscape; 
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ii) next in importance is the primary homoclinic orbit 
to the period 2 orbit in i), (1/3,1/3), which goes 
to (2/3, 1/6), quickly gets into the region of the pe- 
riod two orbit and is difficult to resolve. The par- 
ity and time-reversal symmetric image points are 
also indicated. It turns out that there is an infinite 
set of periodic orbits which approximate this orbit 
more and more closely. Its effects may be present 
simultaneously, and indistinguishable from the ho- 
moclinic orbit itself ]l2|j . The relevant family (set) 
is denoted by (001) O i which is based on a complete 
binary coding of the orbits [^5|. For example, the 
first few periodic orbits of the family are associ- 
ated with the binary codes (00101), (0010101), and 
(001010101). They are also shown in Fig. (3b), in- 
cluding the symmetric image points. In the 1-d 
slice of Fig. (2a) we see this orbit as well; 



principal homoclinic excursion; 
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FIG. 3. Comparison of classical structures in the correla- 
tion at N — 100. Details in the text. 



iii) there is an infinite number of orbits homoclinic to 
the period-2 orbit. They become increasingly more 
complicated. The next associated periodic orbit 
family, (0011)oi, is shown in Fig. (3c), including the 
symmetric image points. This family was noted by 
Saraceno to scar eigenf unctions JlJ] . Also shown in 
this figure is the period-4 along the diagonal lines: 
(3/5,3/5) -» (1/5,4/5) - (2/5,2/5) -> (4/5,1/5). 
Figure (4) shows sections at p a = 3/5,4/5,2/5 to 
highlight this orbit. In Figs. (4a) p a = 3/5 and 
has a local minimum at q a = 3/5; (b) p a — 4/5 
and has a local minimum at q a = 1/5; and (c) 
p a = 2/5 and has a local minimum at q a = 2/5. 
These are marked, to indicate location along q a by 
filled circles. The other minima are due to compet- 
ing nearby structures of the period-2 orbit and its 
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FIG. 4. Sections of the correlation (N = 100) to highlight 
the period-4 orbit, (a) p a — 3/5 and has a local minimum at 
q a = 3/5 (b) p a = 4/5 and has a local minimum at q a = 1/5 
(c) p a — 2/5 and has a local minimum at q a = 2/5. 

iv) the orbit homoclinic to the period-4 orbit included 
in Fig. (3c) with the initial condition (1/5,2/5) 
(and its symmetric partners) is shown in Fig. (3d); 
and 

v) points, such as (0, 1/4), which are homoclinic to the 
fixed points (0,0), (1, 1) also show prominently. 




FIG. 5. The correlation at (a) N = 128 and (b) N = 200. 
Note the sharp features in (a), where the peak height is about 
twice as large as that in (b). 

That these structures are in a sense invariant, i.e. not 
specific to N = 100 is shown in Fig. (5a,b) where the cor- 
relation (absolute value) is shown for N — 128 and 200 
respectively. The phase-space resolution of the correla- 
tion is increasing with N, while the overall magnitude is 
decreasing as 1/N. The peculiar properties of the quan- 
tum bakers map for N equaling a power of two po] , p^| , p9[ 
is tested by N — 128. Here the correlation is "cleaner" 



G 



and the stable and unstable manifolds at 1/4, 1/2, and 
3/4 of the fixed points are clearly visible. The peaks are 
well enunciated as well. Both Figs. (5a, b) have contours 
up to 2/3 peak height, so a direct comparison is mean- 
ingful. Higher N values show more clearly the secondary 
homoclinic orbit to the period-2 orbit. 

We may compare these structures with the inverse par- 
ticipation ratio defined as: 



N-l 



P(a) = IW* 



(23) 



It is illustrated in Fig. (6). It shows marked enhance- 
ments at the period-2 and period-4 (along the symmetry 
lines) orbits, and closer examination reveals all orbits up 
to period-4 are present and one orbit of period-6 along 
the symmetry lines (the diagonals); see Ref. Q for a 
more detailed discussion. 




FIG. 6. The participation ratio for the bakers map. Clas- 
sical structures are present in this quantity as well. 



E. General operators and selective enhancements 

The results so far have dealt with the special case 
A(q) = cos(2nq). It seems natural to suspect that the 
structures highlighted in the correlation are dependent 
on the choice of the operator. This turns out to be true, 
and we show here how this works in the bakers map. We 
reemphasize though that were the eigenstates behaving 
crgodically, the correlations would have been consistent 
with zero to within statistical uncertainties independent 
of the choice of the operator. In this sense, a complete 
view of the extent to which the eigenstates manifest phase 
space localization properties comes only from consider- 
ing both the full phase plane of wave packets and enough 
operators to span roughly the space of possible pertur- 
bations of the energy surface. The flexibility of operator 



choice does provide a means to enhance selectively partic- 
ular features of interest supposing one had a specific lo- 
calization question in mind. As an illustration, note that 
localization about the period-3 orbit barely appeared in 
the contour plot of Fig. (3), and yet, we show below that 
it can be made to show up prominently with other oper- 
ators. 



r 



vy 



FIG. 7. The correlation for the series of operators 
A = cos(2n2 m q). Shown are m = (large dot), m = 1 
(small dot), m — 2 (dashed line). 

Since the case A(q) = sin(27mg) has vanishing cor- 
relations for any integer n due to symmetry, the other 
cases of interest are the higher harmonics of the cosine. 
Therefore consider: 



A(q) — cos(2imq). 



(24) 



If n = 2 m for some positive integer m, a rather remark- 
able scaling property of the quantum bakers map is re- 
vealed that is actually implicit in the way the bakers map 
was originally quantized in |25| . Semiclassically, the cor- 
relations are identical to the case m = 0. For example, 
consider A(q) = cos(47rq). Then the one-step back clas- 
sical correlation becomes identical to the zero-th order 
correlation corresponding to A(q) = cos(27T(7). The corre- 
lations all shift by m in the sense that C(l) — ► C(l + m). 
Thus, there is a kind of scale invariancc in the correlation 
like classical fractals, although this is not self-similarity 
in the same curve. Quantum calculations reflect this in- 
variance to a remarkable degree as seen in Fig. (7) where 
the N = 200 and p a = 1/3 case is shown. 

Other harmonics do weight differently the same local- 
ization effects (classical structures). In Fig. (8), N andp Q 
are taken the same as in Fig. (7). The cases n = 1, 3, 5 are 
all very different from each other, but note that the case 
n = 6 almost coincides with n — 3 for the same reason 
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that powers of two harmonics are nearly same. Thus only 
operators of odd harmonics give the possibility of provid- 
ing new or unique information about the nonergodicity in 
the eigenstates. The period-2 orbit localization is accen- 
tuated at n = 3, since for n — 3m where m is a positive 
integer, cos(67rmg) has a maximum of +1 at q = 2/3, 
whereas for all other integers n, cos(27m<7) = —1/2 at 
the same point. In short the perturbation (or measure- 
ment) is more significant at the location of the period-2 
orbit for n = 3. On the other hand, the case n = 5 
is similar to the fundamental harmonic case at q = 2/3 
where the perturbation is also equal to —1/2. 
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FIG. 8. Correlation for different harmonics A — cos(2nnq) . 
Shown are n — 1 (solid large dots), n = 3 (solid line), n — 5 
(dotted line), and n — 6 (dashed line). 

The case n — 7 is interesting as cos(147rq) has a maxi- 
mum at q = 1/7 which coincides with a period-3 orbit at 
(1/7,4/7). In Fig. (9), we see the correlation (N = 100) 
corresponding to this operator and the dominant struc- 
ture is this period-3 orbit and its symmetric partner. 
Also visible are the stable and unstable manifolds of these 
orbits. In fact, it is the multiples of the 2 m — 1 harmonics 
which selectively highlights the period m orbits. 

Summarizing then, the correlations reflect that the 
bakers map eigenstates are not ergodic, and manifest 
strongly phase space localization properties. There do 
not exist transport barriers such as cantori or diffusive 
dynamics in the bakers map, so whatever localization 
that exists should be due to scarring by the short pe- 
riodic orbits. This is confirmed in the examples shown 
with connections to their homoclinic orbits illustrated as 
well. The perturbation or observable determines the re- 
gions of phase space that will light up in the correlation 
measure. A semiclassical theory predicts reasonably well 
many of these structures. The correlation is semiclassi- 



cally written as a sum of classical correlations that are 
super-exponentially cut off after about half the log-time 
scale. 
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FIG. 9. The correlation for the operator A — cos( 147rg) 
and case N = 100. The highlighted areas are in the region of 
classical period-3 orbits. 



III. THE STANDARD MAP 
A. The map and the mixed phase space regime 

The standard map (a review is found in fl35fl) has many 
complications that can arise in more generic models and 
we turn to their study. It is also an area-preserving, two- 
dimensional map of the cylinder onto itself that may be 
wrapped on a torus. We will consider identical settings of 
the phase space and Hilbert space as for the bakers map 
discussed above. The standard map has a parameter that 
controls the degree of chaos and thus we can study the 
effect of regular regions in phase space, i.e. the generic 
case of mixed dynamics. 

The classical standard map is given by the recursion 

Qi+i = (qi + Pi+i) mod (1) (25) 
Pi+i = (Pi ~ (fc/27r)sin(27r%))mod(l), 

where i is the discrete time. The parameter k is of prin- 
cipal interest and it controls the degree of chaos in the 
map. Classically speaking, an almost complete transi- 
tion to ergodicity and mixing is attained above values 
of k w 5, while the last rotational KAM torus breaks 
around k « .971. 

The quantum map in the discrete position basis is 
given by ||| 
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{n\U\n') = ^=exp(in(n-n') 2 /N) 

x exp (i^- cos(27r(n + a)/N)j . (26) 

The parameter to be varied will be the "kicking strength" 
k, while the phase a = 1/2 for maximal quantum sym- 
metries, and n, n' — 0, . . . , N — 1. 

We use the unitary operator and evaluate the corre- 
lation as in Eq. (J3J) with A(q) — cos(2-7r<7) here as well. 
This corresponds exactly to the level velocity induced by 
a change in the parameter k. In Fig. (10) is shown the ab- 
solute value of the correlation for various values of param- 
eter k. Case (a) corresponds to k — .1 x (2tt) and is dom- 
inated by the KAM curves as the perturbation has not 
yet led to significant chaos. Highlighted is the fixed point 
resonance region at the origin that is initially stable. An 
unstable point is located at the point (1/2,0 or 1). The 
separatrix or the stable and the unstable manifolds of 
this point are aligned along the local ridges seen in the 
correlation. Also the period-2 resonance region is visi- 
ble. Higher resolution not shown here, corresponding to 
higher values of N reveal weakly the period-3 resonance 
as well. Case (b) corresponds to k — .3 x (2ir) while 
(c) and (d) to k = .9 x (2tt) and k = 2.3 x (2tt) re- 
spectively. We note the gradual destruction of the KAM 
tori and the emergence of structures that are dominated 
by hyperbolic orbits. A more detailed classical-quantum 
correspondence is however not attempted here. 





FIG. 10. The quantum correlation for the standard map 
(N = 100), (a) fc/(2ir) = .1, (b) k/(2n) = .3, (c) fc/(2vr) = .9, 
(d) fe/(27r) = 2.3. 

These contour plots do not reveal the difference in the 
magnitude between the correlations in the stable and un- 
stable regions. In Fig. (11), we have plotted the corre- 



lation at the origin (0,0), which is also a fixed point, 
as a function of the parameter. The value fc/(27r) ~ 1 
corresponds to a transition to complete classical chaos 
and is reflected in this plot as erratic and small oscilla- 
tions. The large correlation in the mixed phase space 
regime arises from the non-ergodic nature of the classi- 
cal dynamics. The classical fixed point loses stability at 
k* /(2-7r) = 4/(27r) ~ .63 and this is roughly the region 
at which the correlation starts to dip away from unity 
toward lower values. 




O 0.4 - 
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k/(2 ir) 

FIG. 11. The quantum correlation for the standard map at 
a — (0, 0) as a function of the parameter k (N — 100 is the 
solid line and N = 200 is the dashed one). The dotted line 
is a classical estimate based on Eq. ( p8[ ) and averaged over 
twenty time steps with N = 100. 

The gross features and principal fi behavior in this 
regime is easy to derive in terms of purely classical cor- 
relations as follows: 



C A (a) = (Tr(|a)(a|i(n))) n 



dqdp[\a)(a\] w [A(n)} w 



(27) 



where [.]w is the Weyl-Wigner transform of the operator 
in the brackets and A(n) is the operator after a time n. 
Without worrying about the toral nature of the phase 
space and the Weyl-Wigner transforms, we treat the 
problem as in a plane. This is justified by the use of local- 
ized, Gaussian wave packets. Otherwise, we could imag- 
ine that the Wigner transform of the projector would 
follow from an infinite series of Gaussian states that is 
equivalent to discretizing the Gaussian. We use a nor- 
malized, "circular" Gaussian of width \fh. The approxi- 
mation comes in when we replace [A(n)]v7 by A(f n (q,p)) 
where the latter is the classical function evaluated at the 
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classically iterated point q n — f n (q 7 p). We expect this 
approximation to be valid in the case of regular dynam- 
ics over a much longer time scale than found with chaotic 
dynamics. To a good approximation, 



C A (a) 



dqdp [ — 
irn 



exp 

ij>~p a ) 2 )] A(r(q,p))) 



(28) 



As intuitively expected, there is no principal h depen- 
dence in the correlation since there is a non-zero clas- 
sical limit. At h — 0, we could replace the Gaussian 
forms by S— functions and would get simply Ca(oc) — 
(A(P(q a ,p a ))) n . 

This however vanishes as the classical system becomes 
more ergodic and is no more capable of predicting the 
correlation. Higher-order corrections are needed. It is in 
this regime that we studied the bakers map and found 
that the correlation has a principal part that scales (al- 
most) as h and classical correlations based on periodic 
orbits predict the localization features that arise out of 
quantum interference. 

We return to Fig. (11) to remark on some of these 
properties. Notice that the simple estimate of Eq. fl2S| ) 
performs very well, even as the phase space is becoming 
increasingly chaotic. It is quite unexpected that the oscil- 
lations after the onset of full mixing (around kj (2%) = 1) 
should follow this estimate. However, after the transition 
to chaos the classical estimate will depend on the times 
over which the averaging is done and as this increases the 
estimate would vanish. 



B. Chaotic regime 

We attempt in some measure a semiclassical theory 
for the correlation in the chaotic regime along the lines 
adopted for the quantum bakers map. Of the two ingre- 
dients in Eq. (^) one of them remains the same, namely 
Eq. (||). However the diagonal elements of the propaga- 
tor in Eq. (0) have to be generalized. 

In a semiclassical expression for the matrix ele- 
ments of the propagator as a homoclinic orbit sum is 
given. Although this was derived with the example of 
the billiard in mind, it can be interpreted as a gener- 
alization of Eq. (|?j) for area-preserving, two-dimensional 
maps. We, however, interpret the sum not as a homo- 
clinic orbit sum, but as a periodic orbit sum. To each ho- 
moclinic orbit there is a neighboring periodic orbit that 
we will use instead. This will form the points around 
which the expansions are carried out and the result is 
identical to that in |p2l. Thus we write 



^e?ty{iS 1 /h-iiTv/2)^B j (29) 



where 



B; 




A 1 2hA Q 
+ 2i(m 2 \ Sp 2 - mi2 Sq 2 + SqSp (m 2 2 



[(tr - 2)(Sq 2 + 5p 2 ) (30) 
mil))]}. 



Bj generalizes the Gaussian form (including the prefac- 
tor) in Eq. ([?]) . Again j labels points along the periodic 
orbit and Sq = q.j — q a , dp = pj — p a are as before de- 
viations from the centroid of the wave packet. The two 
dimensional matrix elements, the elements of the 

stability matrix at the periodic point j along the periodic 
orbit 7. The deviations Sq and Sp after n iterations of 
the map are given by: 



Sp n 
Sq n 



run 
m 2 \ 



Sp 
Sq 



(31) 



The invariant is the trace of this matrix that is denoted 
tr. While A = m n + m 22 + *( m 2i — m^), v is a phase 
that will not play a crucial role below. In the case of the 
bakers map rti\i = 77121 = and raw — 2~ n , 77122 = 2" 
uniformly at all points in phase space, as well as v = 0. 
On substitution of this in Eq. (^9 ) we get Eq. (fj]) . 

The dependence on individual matrix elements of the 
stability matrix complicates the use of this formula in 
general. However we note that the Gaussian is effec- 
tively cutting off periodic points that are not close to a 
and therefore we may take the my elements to be the 
stability matrix at this point. In the chaotic regime each 
of the matrix elements grow exponentially with time n. 
Thus we have that exp(— An) my — > constant, where A is 
the Lyapunov exponent. We call this saturated constant 
7n^ j as well. Below we will assume that the exponential 
growth has been factored out of these elements. Also we 
use exp(— An) Aq — > an. The terms inside the exponen- 
tial function in Eq. (|30|) saturate in time n while the 
prefactor goes as exp(— An/2). It follows then that 



B, 



— exp{-\n/2) F{q^ P j) 



where F(qj,pj) is 



exp 



[ 2?Lan 



[(Sq 2 +Sp 2 ) 



(32) 



(33) 



2i(m 2 i Sp 2 - 77112 Sq 2 + SqSp (m 2 2 - mn))] } 



Here the my elements already have the exponential be- 
havior factored out. For example in the case of the bakers 
map 77122 = 1 while all the other elements are zero and 
this gives consistently the approximated Gaussian form 
in Eq. (|l^) . Further steps are identical to the case of the 
bakers map and leads to the generalization of Eq. (|ll"|): 



M 



T l=—M 



(34) 



where the classical correlations are calculated as in Eq. 
(11) with the function F being that in Eq. @. We 
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may then expect all the principal conclusions from the 
study of the bakers map to be carried over, principally 
the decrease in the correlation as %, the correlations being 
cut off after half the log-time scale, and the effects of 
classical orbits. 

More detailed analysis in the lines of the special case 
discussed in case of the bakers map will run into the fol- 
lowing difficulties. First, the rriij elements will depend on 
a in general. Exceptions are uniformly hyperbolic sys- 
tems such as the cat or sawtooth maps (and, of course, 
the bakers map). A second difficulty is that the correla- 
tions have to be evaluated to half the log-time while clas- 
sically iterating the map (analytically) over such times is 
often not possible. The classical correlations that arise 
in the study of rms values of level velocities [^H involved 
correlations that exponentially decreased in time while 
here we are likely to get generalizations of forms such as 
in Eq. ( |l6| ) that will require us to go up to log-times. 
We have calculated the correlations for times 1,-1, and 
2 but will not display them as they are by themselves 
not very useful. A third problem with this form of the 
generalization is that it is not explicitly real. 




12 3 4 

k/(2 it) 

FIG. 12. The quantum correlation for the standard map 
(TV = 100) at a — (0, 0) as a function of the parameter k 
(dashed line). The solid line is the semiclassical estimate. 

We have used Eq. (^4|) and for the my used either those 
calculated at one point in phase space (such as the origin) 
or in fact assumed those that are relevant for the bak- 
ers map. While fine structures are not reproduced, the 
general features are captured equally well in both these 
approaches. To illustrate the quality of the approxima- 
tion we again look at the correlation at the point (0, 0) 
as a function of k in Fig. (12) (as in the previous figure). 
The solid line is the semiclassical prediction based upon 
using the same rriij values at all values of k. It is seen 



that even with these (over) simplifications the semiclassi- 
cal expressions capture much of the oscillations with the 
parameter and the magnitude. 

IV. SUMMARY AND CONCLUSIONS 

We have studied the details of phase space localiza- 
tion present in the quantum time evolution of operators. 
This was related to a measure of localization involving 
the correlation between the level velocities and wavefunc- 
tion intensities. While individual quantum states show 
well known interesting scars of classical orbits, groups of 
states weighted appropriately provide both a convenient 
and important quantity to study scmiclassically. We were 
interested principally in those features whose origins were 
quantum mechanical. The quantities studied had both a 
vanishing classical limit as well as vanishing RMT aver- 
ages. 

We studied simple maps as a way to understand the 
general features that will appear. We found that the 
operator dictated to a large extent which parts of phase 
space will display prominent localization features and fur- 
ther that these localization features are often related to 
classical periodic orbits and their homoclinic structures. 
The time average of the operator for wave packets was 
explicitly related semiclassically to classical correlations. 
These were shown to be cut off super-exponentially after 
half the log-time in the quantum bakers map. Thus the 
localization features in quantum systems associated with 
scars were reproduced using long (periodic) orbits but 
short time correlations. The localization would disappear 
in the classical limit as the magnitude of the quantum 
correlations or time averages are proportional to (scaled) 
H. 

General systems were approached using the quantum 
standard map and complications that would arise were 
discussed. Also the case of mixed phase space was seen 
to be well reproduced by a simple classical argument. 
The generalization to Hamiltonian systems Bjj ] contains 
many of the features and structures are also (not surpris- 
ingly) present in this case. 
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